Mon. Not. R. Astron. Soc. 000, [Tlfl3l C20091 Printed 13 January 2013 



(MN MfeK style file v2.2) 



Rest-frame ultraviolet-to-optical spectral characteristics of 
extremely metal-poor and metal-free galaxies 



o 

< 



o 

U 

6 



(N 
> 
O 

in 
in 

o 



is 



Akio K. Inoue 1 * 

1 College of General Education, Osaka Sangyo University, 3-1-1, Nakagaito, Daito, Osaka 574-8530, Japan 



ABSTRACT 

Finding the first generation of galaxies in the early Universe is the greatest step forward 
for understanding galaxy formation and evolution. For strategic survey of such galaxies 
and interpretation of the obtained data, this paper presents an ultraviolet-to-optical 
spectral model of galaxies with a great care of the nebular emission. In particular, we 
present a machine-readable table of intensities of 119 nebular emission lines from Lya 
to the rest-frame 1 /zm as a function of metallicity from zero to the Solar one. Based 
on the spectral model, we present criteria of equivalent widths of Lya, He n A 1640, 
Ha, H/3, [O III] A5007 to select extremely metal-poor and metal-free galaxies although 
these criteria have uncertainty caused by the Lyman continuum escape fraction and 
the star formation duration. We also present criteria of broad-band colours which will 
be useful to select candidates for spectroscopic follow-up from drop-out galaxies. We 
propose the line intensity ratio of [O in] A5007 to H/3 < 0.1 as the most robust criterion 
for < 1/1000 of the Solar metallicity. This ratio of a galaxy with a few Mq yr _1 at 
z ~ 8 is detectable by spectroscopy with the James Webb Space Telescope within a 
reasonable exposure time. 
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1 INTRODUCTION 

Understanding galaxy formation and evolution is one of the 
most important issues in the modern astronomy. Finding 
the first generation of galaxies is the largest step to solve 
the question, while it is also an open question what the first 
generation of galaxies is. Galaxies in the very early Universe 
should be metal-poor or may be even metal-free. Therefore, 
the most metal-poor galaxies would be the first generation. 

The galaxy with the lowest known metallicity is I Zw 18, 
a blue compact dwarf galaxy in the local Universe. The mea- 
sured gas metallicity (t o be precisely oxyg en abundance) is 
about 1/50 of the Sun jlzotov et aLlll997Tl. If we adopt the 
classi cal Solar metallicity Zq = 0.02 ( Anders fe Grevessd 
1 19891 ). it corresponds to Z = 0.0004. At z ~ 2-3, Lyman 
break galaxies (LBGs) selected by the so-called drop-out 
technique have Z ~ 0.002-0.01 (= 1/10-1/2 Zp,) '(e.g., 
Pettini et all l200ll : lErb et all l200rj : iMannucci etafl 120091 : 



measured metallicity is that in the ionized gas not in the 
stellar atmosphere. 

There are hundreds of stars whose metallicity in the 
atmosphere is found to be extremely low Z < 1 x 



1Q- 5 (= 1/2000 Zq) i n the halo of the Galaxy (e.g., 
iBeers fe Christliebl 120051 ). The mass of such extremely 
metal-poor (EMP) stars in the present epoch is small (~ 1 
Mb), while their high- ma ss counter part should exist in the 
past (|Komiva et al.l 120071 ). Indeed, there are some possible 
signatures suggesting that LBGs and LAEs at z > 3 con- 
tain massive EMP stars or even metal-free stars, the so- 
called Population III (Pop III) stars, of a non- negligible 
fraction in their stellar mass dMalhotra fe Rhoadsl 120021: 
Jimenez fe Haiman 2006: iBouwens et al.ll2010al : llnoue et al.l 



201 lh . 



Erb et al.1 l201Ch . The metallicity measurements for galax- 



ies selected by the strong Lya emission line, Lya emitters 
(LAEs), are still very rare beca use of the difficulty of sp ec- 
troscopy for these faint galaxies. iFinkelstein et all (|201ll ) re- 
port Z < 0.004 (= 1/5 Zq) at z ~ 2. These measurements 
are made by some strong nebular emission lines, so that the 

* E-mail: akinoue@las.osaka-sandai.ac.jp 



To find the first generation of galaxies, we are push- 
ing out the redshift frontier. The current record of redshift 
measured by spectr oscopy is z = 8.55 of UDFy-38135539 
ijLehnert et al.ll2010l ). Like this object, the brand-new Hub- 
ble Space Telescope/ Wide Field Camera 3 (HST/WFC3) 
imaging enabled us to select LBG s at z > 7 a nd even 
at z ~ 10 (e.g., IBouwens etafl l2010bl . I20T1I ). These 
highest-z galaxies may contain more EMP or Pop III 
stars than lower-z LBGs. The very blue ultraviolet (UV) 
colours of the HST/W FC3 LB Gs suggest such a p ossibility 
jBouwens et al.ll2010al ) (but see lDunlop et alj|201lh . Future 
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telescopes such as 30-m class extremely large telescopes and 
James Webb Space Telescope (JWST) will further push out 
the redshift frontier by imaging observations. 

Spectroscopy is finally required to measure metallicity 
of galaxies as well as the precise redshift of them. However, 
it is time-consuming because of the faintness of the target 
galaxies. Thus, any pre-selections for spectroscopy will be 
useful. We will start from a sample selected by the standard 
drop-out technique first. Then, it may be useful if we can 
select EMP or even metal-free candidates only with imaging 
data, in particular, broad-band photometric colours. This 
paper presents such a method. 

In this paper, we will present a spectral model of EMP 
or even metal-free galaxies with a great care of the nebu- 
lar emission (both lines and continuum). Young starburst 
galaxies emit strong Lyman continuum (i.e. hydrogen ion- 
izing continuum with wavelength less than 912 A; here- 



Table 1. Nebular parameters explored in this paper. 
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2009, 2010 
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2010;lRaiter. Schaerer. & Fosburv 


I2010T). This model will be 



useful to discuss the physical nature of very high-z LBGs 
with broad-band data (i.e. so-called spectral energy distri- 
bution [SED] fit) as well as to select EMP or metal-free 
galaxies for follow-up spectroscopy. 

In section 2, we will describe the modelling of nebular 
emission and SED of galaxies. In section 3, we will present 
the resultant spectra from UV to optical in the rest-frame for 
various metallicities and present equivalent widths of emis- 
sion lines as a function of metallicity and star formation 
duration. In section 4, we compare our model with the ob- 
served broad-band colours of z ~ 7-8 galaxies and discuss 
how we find EMP and Pop III galaxies in the future. In the 
final section, we will present a summary of this paper. 

This paper follows a standard ACDM cosmology with 
Qm ~ 0.3, Aa = 0.7, and h — 0.7. All the magnitude is 
described in the AB system. 



2 MODEL 

2.1 Stellar spectra 

The SED of pure stellar populations depends on metallicity 
Z, initial mass function (IMF), star fo rmation history , and 
age. Here, we assume a Salpeter IMF (|Salpeter| [l955) with 
1-100 Mq and a constant star formation with a duration 
of 1, 10, 100, or 500 Myr. The metallicities considered are 
Z = 0.02 (= Z Q ), 0.008, 0.004, 0.0004, 1 x 10~ 5 (EMP), 
1 x 10~ 7 , and (Pop III). SEDs of the former 4 cases are 
generated b y the population synt hesis code starburst99 
version 5.1 (|Leitherer et al.l Il999l). Those for the latter 3 
cases are taken from ISchaererl (|2002l . l2003h . 



2.2 Nebular emission 

2.2.1 Lines 

Emission lines of hydrogen, helium, and some other ma- 
jor elements are taken into account. Based on a large 
grid of photo-ioniz ation models by using CLOUDY 08.00 
jFerland et al.l[T9 985), we have obtained emissivities of emis- 
sion lines relative to H/3. Table 1 shows the parameter space 



Parameter 



Values 



log 10 (Z/Z @ ) 
log 10 U 

lo glo( n H/cm- 3 ) 



-5.3, -3.3, - 
-3.0, -2.5, 
0.0, 0.5, 



-1.7, -0.7, -0.4, 0.0 
-2.0, -1.5, -1.0 
1.0, 1.5, 2.0 



Table 2. Average electron temperatures in the nebulae for vari- 
ous metallicities. 
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of nebulae explored by the CLOUDY calculations in this pa- 
per: 7 metallicities Z, 5 ionization parameters U, and 5 
hydrogen number densities nn. The set of Z is the same 
as that of the stellar spectra adopted (§2.1). For each case 
of Z, we input the corresponding stellar spectrum into the 
code. Since there are 4 cases of the duration of star for- 
mation for each Z, the total number of the model calcu- 
lated is 700 (= 4 x 7 x 5 x 5). The considered U and n H 
are typical ranges of Galactic and ext ra-galactic H n re- 
gions (e.g., lOsterbrock fc Ferlandl fcoOG) and the validity of 
the values will be confirmed by a comparison with obser- 
vations for some strong lines later (Fig. 1). We assume the 
plane-parallel geometry for all the calculations and assume 
the abundance of elements to be the Solar one as the de- 
fault in the code. We consider two cases with and without 
dust in ionized nebulae (see §2.2.3 for details). The electron 
temperatures in the nebulae are calculated in the code with 
energy balance between heating and cooling rates. Table 2 
shows average temperatures for each Z cases which are used 
to calculate H/3 luminosity by equations (1) and (2) later 
and nebular continuum emission in §2.2.2. The higher tem- 
peratures in dusty cases are probably due to photoelectric 
heating by dust. 

We select 119 emission lines between Lya and 1 urn 
in the rest-frame: H Lya, Balmer and Paschen series, and 
He and other elements' lines which have an emissivity more 
than about 1% of H/3 in at least one case calculated (except 
for H lines, for which we adopt all the lines output from the 
CLOUDY: up to n — 24). The list of the lines is presented in 
Appendix. 

For each stellar spectrum (i.e. each Z and star forma- 
tion duration), we calculated 25 set of (U , nn). Among the 
25 cases, line emissivities relative to H/3 vary within about 
an order of magnitude for metal lines but the standard devi- 
ation relative to the average is about 5-25%. The variations 
of H and He lines are much smaller and the relative standard 
deviations are typically < 5%. For calculations in later sec- 
tions, we adopt average line emissivities obtained for each 
set of Z and star formation duration. We have also found 
that the line emissivities are very stable against the change 
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Figure 1. Emissivities relative to H/3 of six strongest metal emission lines as a function of metallicity. The circles are results calculated 
with the code CLOUDY 08.00 (Ferland 1998) in this work: no dust model (filled) and dusty model (open). The squares are empirical 
results by Anders & Fritze-v. Alvensleben (2003). The triangles are the observations of I Zw 18 (northwest and southeast components) 
measured by Izotov et al. (1999). The diamonds are observations of a z = 2.3 galaxy measured by Erb et al. (2010). The dashed lines 
are empirical results by Nagao et al. (2006) and updated by Maiolino et al. (2008). 



of star formation duration if it is larger than 10 Myr. This is 
probably because the saturation of the LyC luminosity and 
spectrum after > 10 Myr star formation. In Appendix, we 
present machine-readable tables of line emissivities relative 
to H/3 as a function of only Z (i.e. averaged over 3 star for- 
mation durations of 10, 100, and 500 Myr) which may be 
useful for future calculations by readers. On the other hand, 
in this paper, we adopt the line emissivities depending each 
star formation duration. 

Figure 1 shows the relative emissivities of 6 strongest 
metal emission lines as a function of metallicity Z (ones 
presented in Appendix, i.e. averaged over the star forma- 
tion durations). The circles are theoretical results of this 
paper by CLOUDY 08.00: no dust case (filled) and dusty 
case (open). The sq uares and dashed lines are empir i cal re - 



lati ons compiled by Anders fc Fritze-v. Alvensleben ([20031 ) 
and iMaiolino et all <|2008h (see also iNagao et al.ll200eh 7 re- 



spectively. The triangles and diamonds are observations of 
a local galaxy, I Zw 18 , (|lzotov et al.| [T999) and a z = 2.3 
galaxy (|Erb et alj|201oh . We find a good overall agreement 
between ours and theirs but an order of magnitude differ- 
ence is found for some cases. This indicates a large uncer- 



tainty of the line emissivities which depend o n the nebular 
physic al parameters such as U and riH- Indeed. [Nagao et al.l 
(2006) suggest Z dependence of U to explain their empir- 
ical relations. Despite such uncertainties and difficulties, 
our emissivity of the stron gest line, [Q ill] A 5007, excel- 
lently agrees with that of INagao et al.l l|2006l ) who com- 
piled the largest sample of galaxie s distributed over th e 
widest range of Z and upda ted by Maiolino et al.l 12 008). 
As p o inted out by literature (IZackrisson. Bergvall. fc Leitej 
l2008l ; ISchaerer fc de Barros! |2009| . l201fi) and also seen in 
later sections, this [O ill] emission line has the largest effect 
on broad-band colours and other metal lines do not affect 
the colours significantly. Therefore, we consider that using 
our theoretical line emissivities (or our choice of U and riu) 
is justified and a large uncertainty of line emissivities does 
not degrade our conclusions based on the [O ill] line. 

To obtain the luminosity of the emission lines based on 
the relative emissivities, we need the luminosity of H/3. We 
assume the following expression: 



7H/3(r c )Q*(l - /esc - /dust) 
Q B (T C ) +Ql(T c ) (/esc + /dust) 



(1) 
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where 7H/3 is the H/3 emission coefficient, Q* is the stellar 
production rate of LyC photons, / csc is the escape fraction of 
the photons, /dust is the fraction of the photons absorbed by 
dust within the ionized gas, qb is the Case B recombination 
rate, and ai is the recombination rate to the ground state. 
The denominator is usually described as just ob- However, 
our expression is correct when the nebular LyC escapes from 
the nebulae and is absorbed by dust with the same proba- 
bility of / esc and /dust as the stellar LyC and the ioni zation 
equilibrium is established in the nebulae l|lnouell2010h . The 
fractions of / eS c and /dust is assumed to be independent of 
wavelength, which is valid for clumpy nebulae as discussed in 
llnouei (|2010h . The LyC escape should be taken into account 
especially for z > 3 because a significant escape wa s de- 
tected at z ~ 3 (|Shaplev et al.ll2006l : llwata et al1l2009l ) and 
such galactic ionizing radiation probably caused the cosmic 
reionization. On the other hand, we omit the LyC absorp- 
tion in this paper (i.e. /dust = 0) because the effect is very 
small for low-metallicity cases, according to a discussion in 
§2.2.3. 

In 7h/3, we have omit the density dependence because we 
only consider densities smaller than 10 2 cm~ 3 which is well 
below a criterion for the sm all density limit (< 10 4 cm~ 3 ; 
lOsterbrock fe F crland 2006?). The electron temperature T c 
dependence of 7h/3 can be approximated to 
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Figure 2. H/3 luminosity per unit star formation rate as a func- 
tion of metallicity. The diamonds, triangles, and squares are the 
cases of the LyC escape fraction / csc = and the constant star 
formation duration of 10, 100, and 500 Myr. The crosses are the 
case of /osc = 0.5 and 100 Myr. The small triangles are the same 
case but obtained by the / osc scaling of equation (3). The arrows 
at the left edge of the panel are the metal-free cases. The LyC 
absorption in nebulae is neglected (i.e. /d us t = 0)- 



7H/3 = 1.23 x 1(T 25 (jjjfyr) erg s" 1 cm 3 , (2) 

which is obtained from Table B.5 in iDopita fc Sutherland! 
|2005t ) for the density of 10 2 cm" 3 . The uncer- 
tainty is less than 4% between 5,000 K and 30,000 
K. This emission coeffic ient is obtained with the 
Case B assumption dDopita fc Sutherland ! 120031 '). Recently 
iRaiter. Schaerer. fc Fosburvl l|201o[ ) have shown that H/3 lu- 
minosity is excellently predicted by the Case B even for ex- 
tremely metal-poor cases. However, they also show that the 
Case B prediction underestimates Lya luminosity. On the 
other hand, our prediction of Lya is based on the ratio of 
Lya to H/3 o btained from the CLOUDY, in whic h the effect 
discussed by IRaiter. Schaerer. fc Fosburvl ||201CJ ) is already 
taken into account. 

Figure 2 shows H/3 luminosity as a function of metallic- 
ity of the stellar population. The luminosity is normalized by 
the star formation rate. The H/3 luminosity increases when 
metallicity decreases because the ionizing photon production 
rate Q* is larger when metallicity is lower. The dependence 
of star formation duration is weak or absent when the dura- 
tion is longer than 10-100 Myr depending on the metallicity. 
The dependence of the escape fraction / esc found in equation 
(1) can be approximated to 



2.2.2 Continuum 

As the nebular continuum, we consider bound-free, free-free, 
and two photon emissions of hydrogen. Any helium con- 
tinuum emission is not taken into account in this paper. 
While this simplification underestimates the nebular con- 
tinuum emission, the helium continuum is not important 
very much because it is weak and negligible in fact (see e.g., 
Ilno'uell20ich . Continua of metal elements are not taken into 
account, either. The luminosity density of hydrogen nebular 
continuum at the frequency v is given by the very similar 
form to equation (1) as 

L _ 7v{T c )Q,(l - / C5C - /dust) 

ces{T c ) + Qi(T e )(/esc + fa ust I 

where 7„ is the emission coefficient of the continuum. 
The volume emissivity y v and the recom bination rate s 
qb and Qi are calculated a s des cribed in llnouei i|2010l ). 
IRaiter. Schaerer. fc Fosburvl (|201dh hav e shown tha t the 
Case B assumption, which is adopted in llnoud (I2OIOI L un- 
derestimates the nebular two-photon emission for the stellar 
effective temperature of > 5 x 10 4 K and the density of < 10 2 
cm -3 . This point causes an underestimation of the nebular 
continuum for A < 0.2 /im, while the effect is negligible for 
A > 0.2 /jm where the nebular continuum is dominated by 
the bound-free emission. 



where the factor 0.6 in the denominator comes from ai/aB 
for T c = 10 4 K. Note that the scaling is not just 1 — / CS c 
as usually assumed when the nebular LyC can escape. This 
scaling is very nice as shown by the comparison of the crosses 
with the triangles in Figure 2. 



2.2.3 Dust 

There are two effects of dust; one is attenuation (or ex- 
tinction) of radiation by dust in the interstellar medium 
(mainly in the outside of ionized nebu lae). This effect is o f- 
ten described with the Calzetti law (|Calzetti et al.l l2000). 
and we do so in §4, although some studi es have suggested 
a different attenuatio n law in high- z (e.g.. lSiana et al.ll2009l ; 
iGallerani et~al]|2010h . 
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Figure 3. Rest-frame model spectra of (a) Pop III (Z = 0) galaxies and (b) moderate sub-solar metallicity (Z = 0.004: log lo (Z/Z0) = 
—0.7) galaxies. In each panel, the lines correspond to 10, 100, and 500 Myr of constant star formation (bottom to top). LyC escape 
fractions / csc =0.5 and are assumed for the solid and the dotted lines, respectively. The emission line width is assumed to be 300 km 
s _1 as an example. No attenuation by dust and IGM is not included. 



The other effect is LyC absorption by dust within ion- 
ized nebulae (i.e. /dust in eqs. [1] and [3]). Dust competes 
with hydrogen to absorb LyC in nebulae. Indeed, a half of 
LyC is directly absorbed by dust in H H regions in the Milky 
Way and some local group galaxies l|lnoue et al.ll200ll ; llnoud 
l200ll ). The amount o f LyC absor ption depends metallicity 
and dust-to-gas ratio (|lnoudl200ll ). To see this effect, we ran 
CLOUDY with the interstellar dust of the code (i.e. Milky Way 
dust). The amount of dust was simply assumed to be scaled 
by metallicity. We also assumed that dust was well mixed 
with gas, while a cen tral dustless cavity is likely to exist in 
nebulae (Inouc 2002). The uniform dust distribution may 
overestimate the effect of LyC absorption. As found from 
Figure 1, the dust effect on emission lines results in larger 
emissivities relative to H,3 in most cases, because LyC ab- 
sorption reduces H/3 emissivity. However, this effect is not 
very large, especially for low-metallicity cases. Therefore, we 
omit LyC absorption (i.e. /dust = 0) and adopt only no dust 
cases in the following sections. 



3 RESULT 

3.1 Rest- frame UV-to-optical spectrum 

Figure 3 shows resultant model spectra in the rest-frame for 
two cases of Z: (a) Z = (Pop III) and (b) Z = 0.004 
(\og 10 (Z / Zq) = —0.7). A constant star formation rate is 
assumed and the duration of star formation is 10, 100, or 
500 Myr (bottom to top in each panel). The LyC escape 
fraction / csc is assumed to be 0.5 (solid) or (dotted) for 
both stellar and nebular continua. We do not consider any 
attenuation by dust and IGM in this figure. Note that the 
vertical axis is normalized by a unit star formation rate. 

In both cases of Z, many H recombination lines such as 
Lycv, Ha, H/3 are highly visible. In addition, some He lines 
such as He n A1640 and He I A5876 are also remarkable in 
the Pop III case. In the Z — 0.004 case, many metal emission 
lines appear. Especially, metal emission lines shown in Fig- 
ure 1 such as [O in] A4959/A5007 lines are very strong. These 
lines affect even the broad-band colours as found laterQ 



1 A spectral dip just longward of the Paschen limit (8204 A) is 
not real but an artifact caused by the lack of Paschen series lines 
higher than n = 24. 



6 A. K. Inoue 




10.0 



0.2 0.4 0.6 O.f 

Rest wavelength [jtini] 



Figure 4. Rest-frame spectra of luminosity density ratio of nebu- 
lar emission to stellar emission: (a) Pop III and (b) moderate sub- 
solar metallicity (Z = 0.004: log 10 (Z/Z(T,) = —0.7). In each panel, 
the solid lines correspond to the LyC escape fraction / esc = 0.5 
and the dotted lines correspond to / csc = 0. For each type of the 
lines, the top to bottom correspond to 10, 100, and 500 Myr of 
constant star formation. The dashed horizontal lines indicate the 
ratio to be unity (i.e. equal contribution from stellar and nebular 
emissions). 



The nebular contribution to the continuum is more sig- 
nificant in the Z = case than in the Z = 0.004 case. The 
Balmer jump in the bound-free continuum is easily recog- 
nised in the Z — and the duration of 10 Myr case. This 
point is shown more clearly in Figure 4 which shows spectra 
of luminosity density ratio of nebular to stellar emissions. 
The lines top to bottom correspond to the duration of 10, 
100, and 500 Myr. As the duration of star formation in- 
creases, the ratio decreases. In addition, we see a stronger 
nebular contribution for a longer wavelength. 

There are three spectral jump in the nebular continuum 
within the wavelength range shown: Lyman, Balmer, and 
Paschen jump of the bound-free continuum. In particular, 
the Balmer jump at 3646 A is strong and can be found in 
the total (stellar+nebular) spectra shown in Figure 3. For 
example, the Pop III case with 10 Myr duration in Figure 3 
shows a factor of ~ 2 Balmer jump. Note that the 'Balmer 
jump' is a sudden decrease of the continuum level towards 
longer wavelength. This is opposite to the 'Balmer break' 
which is a sudden increase of the continuum towards longer 
wavelength found in spectra of older stellar populations. For 
example, the 100 and 500 Myr cases in Figure 3 (b) show a 
prominent Balmer break. 

The strength of the Balmer jump depends on the dura- 
tion of star formation and metallicity. This point is shown 
in Figure 5 more clearly: the nebular-to-stellar ratio at just 
shortward of the Balmer limit as a function of the dura- 
tion. The LyC escape fraction of / csc = is assumed in 
Figure 5, and thus, the ratio shown is the maximum. The 
/eac dependence of the nebular-to-stellar ratio is simply the 
same as equation (3). For all the 7 cases of metallicity, the 
nebular contribution, or the strength of the Balmer jump, 
monotonically decreases as the duration increases. The neb- 
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Figure 5. Luminosity density ratio of nebular to stellar emis- 
sions at 3640 A, very close to the Balmer limit (3646 A), as a 
function of the duration of star formation. There are 7 metallic- 
ity cases: log 1(j (Z/ZQ) = oo (plus), —5.3 (asterisk), —3.3 (circle), 
— 1.7 (diamond), —0.7 (triangle), —0.4 (square), and 0.0 (cross). 
No escape of the LyC is assumed (/ eS c = 0). 



ular contribution also decreases for higher metallicity. For 
> 10 Myr duration, only EMP (log 10 (Z/Z Q ) < -3.3) or 
metal-free cases exceed the ratio of unity, i.e. stronger nebu- 
lar emission. Thus, a prominent Balmer jump appears only 
for these cases, while normally sub-solar metallicity (i.e. 
log 10 (Z/Z(s) —1-7) cases expect a weaker Balmer jump. 
Therefore, the strong Balmer jump can be a signature of 
EMP or Pop III stellar populations. 

In Figure 4, we can find another strong jump at the Ly- 
man limit. This feature is also caused by H bound-free emis- 
sion a nd appears only when / csc > 0. As proposed by llnouel 
(2010), this Lyman jump (or 'bump') can be also useful as a 
signature of EMP or Pop III stellar populations because the 
j ump becom es strong enough onl y when EMP or metal-free 
l|lnouell2010l ). Ilnoue et all ll201lfl show ed that strong LyC 
emission detected bv llwata et alj |2009) from z ~ 3.1 LAEs 
can be attributed to the Lyman jump by EMP or Pop III 
stars. However, thi s feature is easily obscured by IGM at- 
tenuation at 2 > 4 (llnoud koiO) and we cannot use it to find 
EMP or Pop III galaxies at very high-z. 



3.2 Equivalent width of emission lines 

The equivalent width (EW) of an emission line is denned as 



EW lin 



-^- / line(/csc) 



(5) 



where Lii no is the line luminosity, L x is the stellar lu- 
minosity density at the line wavelength \u ae , and is 
the nebular continuum luminosity density at Aii ne - As ex- 
plicitly expressed in the equation, Ln no and L^r depend 
on the LyC escape fraction / esc . The dependence of these 
two terms is the same as equation (3). Thus, the / csc depen- 
dence of EW has two extreme cases: (a) the same form as 
equation (3) when L* x ^> L^ b and (b) independent of 



/ csc when L* x < L" 



Therefore, the / OS c dependence of 



EW is not very simple when the nebular continuum is rel- 
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Figure 6. Rest-frame equivalent width of Lyo emission line as 
a function of the duration of star formation. The symbols corre- 
spond to 7 metallicities same as Fig. 5. No escape of the LyC is 
assumed (/ CS c = 0). 



atively strong. Such cases happen at the rest-frame optical 
when the star formation duration is very short (~ 1 Myr) 
or when metallicity is extremely low (Z < 10 -5 ) as shown 
in the previous subsection. 



3.2.1 Lya 

Figure 6 shows the rest-frame EW of Lya as a function 
of the duration of star formation. The LyC escape fraction 
/esc = is assumed. We can scale the EW for other / csc by 
the /osc dependence same as equation (3) excellently well 
because the nebular continuum at 1216 A is weak enough. 
We see that the EW decreases monotonically as the dura- 
tion becomes longer. We also see that the EW is smaller as 
met allicit y is higher. These are already reported in litera- 
ture (e.g. , ISchaererl l2002 . 200 j; iRaiter. Schaerer. fc Fosburvl 
2010). Note that our EWs are quantitatively very consis- 
te nt with those with the same IMF a nd a similar duration 
in lRaiter. Schaerer. fc Fosburvl (|2010l l. 

We expect the maximum EW for galaxies with Z ^ 
0.0004 (log w (Z/Z Q ) = -1.7) to be 460 A, or to be 230 A if 
we consider only the age larger than 10 Myr (Note that the 
probability to observe a galaxy as young as 1 Myr is gen- 
erally small because of its short time). Therefore, we may 
conclude that galaxies with an EW > 230 A (or > 460 
A for more conservative ) are EMP or Pop III (see also 
iMalhotra fc Rh oads 2002]). However, Lya transfer in the in- 
terstellar medium ( ISM) is_complex, and sometimes, it en- 
hances the EW (Neufcld 1991J). Thus, a galaxy with higher 
metallicity may have a boosted EW which exceeds the crite- 
rion, and then, may be identified as an EMP/Pop III galaxy. 
On the other hand, it is worth noting that the EW becomes 
smaller than the criterion even for EMP/Pop III galaxies 
if /csc > 0. For example, we find that the EW is 65 A for 
Z = and 10 Myr constant star formation when / csc = 0.9. 
Therefore, Lya EW is intriguing but we need other signa- 
tures simultaneously to conclude a galaxy to be EMP/Pop 
III. 




Age [Myr] 

Figure 7. Same as Fig. 6, but for He II A1640. 
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Figure 8. Same as Fig. 6, but for Ha. 

3.2.2 He n 

Figure 7 shows the rest-frame EW of He II A1640. We see 
that the EW > 1 A is realized only when metallicity is zero 
(i.e. Pop III). Therefore, this em ission line is proposed as the 
signature of Pop III stars (e.g.. ISchaerer]|2002j . 120031 '). How- 
ever, this line may be too weak to be detected, and what is 
worse, the EW becomes even smaller than those in Figure 7 
if /csc > 0. The scaling of / CS c in equation (3) is good for 
> 10 Myr cases because the nebular continuum contribu- 
tion is not very large at the wavelength of the line for the 
cases. Even if we detect the He 11 line, there is another issue 
that the line becomes observable from Wolf-Rayet stars with 
normal metallicity. In fact, the line has been detected from 
z = 2-3 LBGs but its origin is attributed to the stell ar winds 
of these stars (|Shaplev et alj|200tj ; lErb et al.ll20ld '). There- 
fore, using the He 11 line needs mo re careful consid erations, 
for example, about the line width l|Erb et al.ll2O10T ). 

3.2.3 Ha and Hf3 

Figures 8 and 9 show the rest-frame EWs of Ha and H/3, 
respectively. These two cases show qualitatively similar re- 
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Figure 9. Same as Fig. 6, but for H/3. 



suits: the EW becomes smaller as metallicity is higher or 
the duration is longer. For ~ 10-100 Myr durations, we can 
separate EMP and Pop III cases from other higher metal- 
licity cases by criteria, for example, EW(Ha) > 1900 A and 
EW(H/3) > 300 A. If we wish to remove very young galax- 
ies with higher metallicities, the criteria should be EW(Hq) 
> 3200 A and EW(H/3) > 540 A. Unlike Lya, Ha and 
H/3 photons do not undergo resonant scattering. Therefore, 
there are any contaminants of galaxies with higher metal- 
licities in the sample selected by these criteria, except for 
AGNs. In this sense, the Ha and H/3 EWs are promising 
tool to find primordial galaxies in the near future. 

Note that we will miss EMP and Pop III galaxies with 
/esc > by these criteria. The / osc dependence of the EWs is 
not simple because the nebular continuum significantly con- 
tributes to the total continuum at Ha and H/3, especially 
for the EMP and Pop III cases. The scaling of equation (3) 
overestimates the EWs by a factor of 2 or more, unfortu- 
nately. We need full calculations in order to obtain the EWs 
for other / csc with higher accuracy. 



3.2.4 [O in] 

Figure 10 shows the rest-frame EW of [O in] A5007 which 
is the strongest metal emission line in our expected spectra. 
As shown in Figure 1, the metallicity dependence of this line 
is not monotonous. The line strength becomes maximum at 
around log 10 (Z/Z@) — —0.7. For higher or lower than the 
metallicity, the EW becomes lower. We expect that EMP 
and Pop III galaxies have the EW < 20 A. If we consider 
/csc > 0, the EW becomes even lower. Thus, all EMP and 
Pop III galaxies satisfy this criterion. On the other hand, the 
EW for a higher metallicity will also satisfy the criterion if 
/csc > 0. This will be contaminant. In addition, very old 
quiescent galaxies with higher metallicities, if they exist at 
high-z, are also contaminant. 
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Figure 10. Same as Fig. 6, but for [O III] A5007. 
4 DISCUSSION 

4.1 Rest-UV colours of high-z galaxies with 
HST/WFC3 

4.1.1 z~7 

iBouwens et al.l |2010af ) have reported very blue UV colours 
of z ~ 7 less luminous LBG s found in the ultra- deep sur- 
vey with HST/WFC3, while iDunlop et all (1201 ll) recently 
challen ged their finding. The argument by IBouwens et al.l 
l|2010ah is that the UV slope /3 = -3 found in the LBGs with 
Muv = —19 to —18 AB indicates extremely low metallicity 
as log 1Q (Z/Z@) ^ —3.3 and large escape fraction of the LyC 
as /esc > 0.3. Let us examine this argument with our model 
spectra. 

Figure 11 shows a comparison o f our model with the ob - 
served UV slope (or J— H colour) bv lBouwens et all (|2010al V 
For each metallicity, we consider four cases of the star for- 
mation duration (1, 10, 100, and 500 Myr) and three cases of 
the LyC escape fraction (0, 0.5, and 1). We find that almost 
all the EMP and Pop III cases with log 10 (Z/Z s ) < -3.3 
reside within a 2-cr range of the observed /3 = —3.0 ± 0.2 
of the less luminous LBGs. The most plausible duration 
and /csc among the EMP and Pop III cases are the du- 
ration of < 100 Myr and / esc = 0.5-1. For the metallic- 
ity log lo (.Z/Z0) ^ —1.7, the observed ft is reproduced only 
when the duration is ~ 1 Myr, while the 10 Myr cases of 
\og 1Q (Z / Zq) = —1.7 and —0.7 still remain within the up- 
per 2-cr range if / CS c > 0.5. It is unlikely to observe galax- 
ies with ~ 1 Myr age because of the short time to observe 
them. Furthermore, the higher metallicity cases are probably 
reddened by dust, which makes more difficult to reproduce 
the observed ft with such highe r metall icities. Therefore, we 
agree with the IBouwens et all (|2010al Vs argument if their 
very blue j3 is real. 

On the other hand, IBouwens et"afl (1201031 ) have re- 
ported the UV slope /3 = -2.0 ± 0.2 for luminous LBGs 
(Muv = -21 to -20 AB) at z ~ 7. Without dust redden- 
ing, we need the Solar metallicity and the star formation 
duration of > 100 Myr to reach the /3. If we apply a small 
amount (e.g., Eb-v = 0.1) of the dust reddening to the 
model, lower metallicity cases can reach the /3. However, the 
EMP and Pop III cases are unlikely to have dust enough 
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Figure 11. HST/WFC3 J125 - #160 colours of z = 7 galaxies 
for seven metallicities. Each metallicity, there are four cases of 
the star formation duration: 1 Myr (circle), 10 Myr (triangle), 
100 Myr (square), and 500 Myr (diamond). The metal- free case 
are indicated by arrows at the left edge of the panel. The four 
models of each metallicity are horizontally shifted each other for 
the display purpose. Each model has three cases of the LyC escape 
fraction / e s C : the upper open symbols for / csc = 0, the middle 
filled symbols for / esc = 0.5, and the lower open symbols for 
/csc = 1. The vertical arrow is the dust reddening vector based 
on Calzetti et al. (2000), which is likely to be applied for higher 
metallicity cases. At the right edge of the panel, the UV slope 
based on the conversion by equation (1) in Bouwens et al. (2010a) 
is shown. The two shaded regions are the observed ranges of 
for luminous and less luminous galaxies at z ~ 7 reported by 
Bouwens et al. (2010a); the thick and thin shades indicate the 
ranges within 1-cr and 2-a uncertainties, respectively. 



0.4 
0.2 



m 
< 



-0.6 
0.8 



n — 1 — 1 — 1 — 1 — 1 — r 



n — 1 — r 




2 = 

Constant SF 

G»G1 Myr _ , ... , 

A*A10 Myr f^ 26 " 1 }?* 

W>1100 Myr K^b-v-V- 1 ) 
O*O500 Myr 

9 



_i I 1 1 1 I 1 1 1 I 1 1 i_ 



-->( 



-6 -4 -2 
logio(^Ao) 



Figure 12. Same as Fig. 11, but for z = 8 galaxies. Note the much 
wider range of the colours than Fig. 11. The crosses with vertical 
error-bars near the right edge of the panel are the eight robust 
sample of z ~ 8 galaxies compiled by Taniguchi et al. (2010). 
Their positions along the horizontal axis are meaningless and just 
for the sake of showing. 



and to reach the observed /?. Therefore, the luminous LBGs 
at z ~ 7 probably have metallicity larger than O.OIZq. 

4.1.2 z~8 



iTanieuchi et all <|2010h have presented a ro bust sample of 
galaxi e s at z ~ 8 detected with HST/WFC 3 (Bouwen s et al 



2010bl; iBunker et all l20ld; iMcLure et al.1 lipid ; lYan et~. 
20ld ; iFinkelstein et alj |2010| ). Let us compare the UV 
colours of the galaxies with our model. Figure 12 shows the 
comparison. However, a large uncertainty and variance on 
the observed colour make it difficult for us to derive any 
implication from the comparison. Thus, we comment just 
one thing about the bluest two objects which are much 
bluer than our bluest model (i.e. EMP and Pop III cases 
with /csc = 1). This may indicate the presence of a strong 
Lyq emission line in the J125 band of the two objects 
(|Taniguchi et al.ll20ld ). Indeed, the line enters into the J125 
band if the object is located at z ^ 8.1. Note that the colour 
at z = 8 shown in Figure 12 does not include the effect of 
the line. 



4.2 Rest-optical colours of high-z galaxies with 
JWST/NIRCAM 

As discussed in section 3.1, there are the Balmer jump and 
some strong emission lines in the rest-frame optical. The 
observed wavelength of these features from high-z becomes 
near- and mid-infrared where ground-based observations are 
difficult. The forthcoming JWST will cover the wavelength 
from the space. Let us examine the effect of the optical fea- 
tures on broad-band colours observed with the Near Infrared 
Camera (NIRCAM) on the JWST. In particular, we focus on 
F244W— F444W colour because the Balmer jump at z ~ 8 
comes between the two bands0 

Figure 13 shows the redshift dependence of the colour. 
Although we show only four combinations of metallicity and 
star formation duration as indicated in the panels, the trend 
of other cases is similar to these cases. We find that the neb- 
ular emission makes the colour redder than the pure stellar 
one. The degree of the reddening increases as the LyC escape 
fraction / esc decreases. This reddening is mainly caused by 
the emission lines like Ha, H/3, and [O ill] as found by the 
comparison between the cases with and without emission 
lines which are indicated by the filled and open circles, re- 
spectively. In fact, the nebular continuum makes the colour 
bluer than the stellar one for z ~ 8 galaxies because of the 
Balmer jump. However, the reddening by strong H/3 and [O 
in] emission lines in the F444W band overcomes the bluing 
by the Balmer jump. We also find that for a longer star for- 
mation duration, the reddening by the nebular emission is 
smaller. In the same time, the stellar colour is redder be- 
cause of the Balmer break of older stars. As a result, the 
colour with the nebular emission (with the both of lines and 
continuum) does not depend on the duration strongly. 

The difference of metallicity appears in the colour. At 
2 ~ 6, Ha is in the F444W band, and then, F277W-F444W 
> for the Pop III cases with / csc > 0.5 but the Z = 0.004 

2 The wavelength coverages of the two broad-bands taken from 
http : / / ircamera.as . arizona.edu /nircam / 
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Figure 13. JWST/NIRCAM F277W-F444W colours of galaxies as a function of redshift: (a) Pop III (Z = 0) galaxies with 10 Myr 
constant star formation, (b) Z = 0.004 (log lo (Z/Z0) = —0.7) galaxies with 10 Myr star formation, (c) Pop III case but 100 Myr star 
formation, and (d) Z = 0.004 case but 100 Myr star formation. Four LyC escape fractions / esc are assumed: / esc = for the filled circles, 
/esc = 0.5 for the triangles, / csc = 0.9 for the squares, and / esc = 1 (i.e. pure stellar spectrum) for the diamonds. The open circles are 
the case with / esc = but no emission lines (i.e. stellar+nebular continua). 



cases are almost always F277W-F444W < 0. At z ~ 8, 
the Z = 0.004 cases show F277W-F444W > due to 
[O in] in the F444W band, but the Pop III cases show 
F277W-F444W < thanks to the lack of the [O in] line. 
Figure 14 shows the metallicity dependence of the colour 
for z ~ 8 galaxies more in detail. As found from the figure, 
the EMP and Pop III cases with log 10 (Z/Z @ ) < -3.3 ex- 
pect F277W-F444W < 0. On the other hand, the higher 
metallicity cases expect F277W-F444W > 0, except for the 
cases of very young (< 10 Myr) and / csc > 0.5. Note that 
the colour of the higher metallicity cases is the lower limit. 
These galaxies probably have dust which reddens the colour 
as indicated by the upper arrow in the panel. Therefore, a 
colour criterion of F277W— F444W < may be useful to 
select EMP and Pop III candidates from galaxies at z ~ 8. 



4.3 Signatures of primordial galaxies 

We discuss what signatures of EMP and Pop III galax- 
ies we can detect in the near future. First, we present the 
selection criteria for EMP/Pop III candidates with broad- 
band colours. Then, we present more robust signatures to be 
searched by follow-up spectroscopy. Table 3 is a summary of 
the criteria. 
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Figure 14. JWST/NIRCAM F277W-F444W colours of z = 8 
galaxies as a function of metallicity. The meaning of the symbols 
are the same as Fig. 11. The horizontal dashed line is a criterion 
to separate the EMP cases from other higher metallicity cases. 
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Table 3. Proposed criteria to find extremely metal-poor and 
metal-free galaxies. 



Broad-band colours 


HST/WFC3 


J125 - #160 < -0.15 for z ~ 7-8 


JWST/NIRCAM 


F277W- 


-F444W > for z ~ 6 




F277W- 


-F444W < for z ~ 8 


Rest-frame EW (A) 


Lyd 




> 230 


Ha 




> 1900 


H/3 




> 300 


He II A1640 




> 1 


[O in] A5007 




< 20 


Metal/hydrogen ratio 


[O m] A5007/H/3 




< 0.1 



4-3.1 Broad-band colours 

As discussed with Figure 11, the rest-frame UV colours can 
be u sed as the signature fo r EMP and Pop III galaxies (see 
also iBouwens et al.ll2010al b However, the colour difference 
between EMP/Pop III and higher metallicity cases is rela- 
tively small, ~ 0.1-0.2 mag. Thus, the observational uncer- 
tainty for individual objects prevents us from distinguishing 
them, e.g., Figure 12. If future facilities suppress the uncer- 
tainty of the colour, the UV colours will be useful as the 
signature for EMP/Pop III stars in individual galaxies. A 
proposed criterion is J125 — Hieo < —0.15 for z ~ 7-8. 

The discussion with Figure 14 suggests the possibility to 
use the rest-frame optical colours as a signature of EMP and 
Pop III galaxies. The colour difference between EMP/Pop 
III and higher metallicities is larger than those of the UV 
colours: ~ 0.2-0.4 mag. In this respect, the optical colours 
is better than the UV colours. However, we need space tele- 
scopes in order to measure the optical colours of very high- 
z galaxies. The JWST and a Japanese project, Wide-field 
Imaging Surveyor for High-redshift (WISHJ3 will be useful. 
Proposed criteria are F277W-F444W > for 2 ~ 6 and 
F277W-F444W < for z ~ 8. 

4-3.2 EW of emission lines 

If emission lines characterising EMP and Pop III galaxies are 
strong enough to be detected, the lines are very useful as the 
signature of these primordial galaxies. The H e 11 A1640 is the 
most discussed feature for Pop III stars fe.g.. lSchaererlr2002l . 
2003). As shown in Figure 7, this line is indeed special fea- 
ture for the metal-free case. However, there are two problems 
to use the line. One is that the strength of the line is not 
very strong. The expected EW is about 1-3 A in the rest- 
frame for the star formation duration longer than 10 Myr. 
This corresponds to the observed EW of < 30 A for z < 9. 
The other is that the strength becomes even weaker if the 
LyC escape fraction / csc > 0. In addition, there is contami- 
nation of t he He 11 line from Wolf-Rayet stars with normal 
metallicity (IShaplev et al.ll2006l ; lErb et al.ll2010r ). Therefore, 
we need other signatures to confirm the galaxies selected by 



the line to be truly Pop III. In any case, we can use the line 
to select the candidates. A proposed criterion is EW(He 11) 
> 1 A in the rest-frame. 

As discussed in section 3.2, EWs of Lya, Ha, and H/3 
can be used as a signature of EMP and Pop III galaxies. 
However, if we adopt a lower limit on the EWs to select the 
primordial galaxies based on the estimations with / esc = 0, 
we will miss some of them which have / esc > 0. In ad- 
dition, Lya photons undergoes resonant scattering in the 
ISM, and someti mes, the EW with higher metallicities is 
boosted by dust (|NeufeIdlll99ir i. and may exceed the crite- 
rion and become contaminant. There is no such case with 
Ha and H/3. Thus, these lines are more reliable signatures. 
Proposed criteria are EW(Lya) > 230 A, EW(Ha) > 1900 
A, and EW(H/3) > 300 A in the rest-frame, or to avoid 
very young (~ 1 Myr) galaxies with higher metallicities, 
EW(Lya) > 460 A, EW(Ha) > 3200 A, and EW(H/3) > 540 
A in the rest-frame. 

With the EW of [O in] A5007, we can select all the pri- 
mordial galaxies by an upper limit on the EW. However, in 
this case, we will have contaminants of galaxies with higher 
metallicities and / GSC > 0. A proposed criterion is EW(0 ill) 
< 20 A in the rest-frame. 

4-3.3 [O III//H/3 line ratio 

Metal-to-hydrogen line ratio is free from the uncertainty of 
/csc, unlike the EWs. As found from Figure 1, the ratio of 
[O ill] A5007/H/3 < 0.1 if log 10 (Z/Z Q ) < -3. Although the 
ratio also becomes < 0.1 if lo g ln (Z/Zp)) > 0.6 ba sed on the 
empirical relation reported by Nag ao et al.l ([2006 ), such very 
high metallicity galaxies would not exist at very high-z, or 
if they exist, itself is very interesting. Therefore, searching 
a galaxy with the ratio < 0.1 is very attractive. 

Based on Figure 2, the H/3 luminosity per a unit star 
formation rate of EMP/Pop III galaxies is - 20-30 x 10 40 
erg s" 1 (Mq yr" 1 )" 1 if the star formation duration is longer 
than 10 Myr and / csc = 0. Thus, the expected upper limit 
on the [O ill] luminosity of EMP/Pop III galaxies is < 2- 
3x 10 40 erg s" 1 (M© yr" 1 )" 1 . This corresponds to the [O ill] 
line flux < 3 x 10" 20 erg s" 1 cm" 2 (M© yr" 1 )" 1 for a z = 8 
source. The expected line sensitivity of the Near Infrared 
Spectrograph (NIRSpec) on the JWST is 3 x 10" 19 erg s" 1 
cm" 2 at 4.5 /um at S/N=10 with R=1000 mode in 100,000 
s exposure time0 Therefore, we can reach [O lll]/H/3 < 0.1 
at 3-cr for a 10 (or 3) Af© yr" 1 galaxy at z = 8 in 10,000 
(100,000) s with the JWST/NIRSpec. 

4-3.4 B aimer jump 

As discussed in previous sections, we expect a significant 
Balmer jump in spectra of EMP/Pop III galaxies. Let us 
examine if we detect this feature with spectroscopy. Based 
on Figure 3 (a), we expect the continuum level around the 
Balmer jump (3646 A in the rest-frame) to be (2-10) x 10 27 
erg s" 1 Hz" 1 (M© yr" 1 )" 1 , depending on / osc and the star 
formation duration. Thus, we expect an observed flux den- 
sity of (2-10) njy (M© yr" 1 )" 1 for z = 8 galaxies. On the 
other hand, the expected sensitivity of the JWST/NIRSpec 



3 http:/ /www. wishmission.org/en/index. html 



4 http: / /www. stsci.edu/jwst/instruments/nirspec/sensitivity/R1000_line.pdf 
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is about 200 njy at 3-4 fj,m at S/N = 10 with R=100 mode 
in 10,000 s exposure timeQ Therefore, we can detect contin- 
uum of a 10 Mq yr _1 galaxy at z = 8 in 100,000 s exposure 
with very high significance of S/N = 3-15. Namely, we can 
detect the Balmer jump of primordial galaxies at z = 8 with 
the JWST/NIRSpec. 



5 SUMMARY 

This paper presents a spectral model from UV to optical 
in the rest-frame of galaxies with various metallicities. A 
special feature of the model is the nebular emission of lines 
and continuum. We have calculated intensities of 119 emis- 
sion lines from H Lya to 1 /im in the rest-f rame with the 
publi c photo-ionization code cloudy 08.00 (|Ferland et alj 
Il998l ). The stellar spec tra input into the co de are gener- 
ated by starburst99 (|Leitherer et al. 199 9J) fo r meta llic- 
ity Z > 1/50Zq or taken from ISchaereil l|2002l . 120031 ) for 
EMP and metal-free cases. We input the stellar spectrum 
with the same metallicity as that in the nebular gas into the 
photo-ionization code. After exploring a wide range of neb- 
ular parameters (Table 1) which are appropriate to real H II 
regions, we have derived average intensities of 119 emission 
lines relative to H/3 (Figure 1 for the 6 strongest metal lines) 
and temperatures of the nebular gas (Table 2) as a function 
of metallicity. The emission line intensities are presented in 
Appendix as a machine-readable form. 

The resultant spectra of galaxies show some interesting 
features in the rest-frame optical (Figure 3). In the EMP and 
metal-free cases, we find a strong Balmer jump, which is a 
spectral plummet towards longer wavelength, by the nebu- 
lar bound-free continuum (Figures 4 and 5). As a result, 
broad-band colours straddling the Balmer jump becomes 
bluer than the stellar one (see Figure 13). However, strong 
Balmer emission lines can fill in the jump and can make 
the colours even redder than the stellar one (Figure 13). In 
higher metallicity cases, some emission lines of oxygen like 
[O ill] are so strong that they redden the broad-band colours 
(Figure 13; see also e.g.. ISchaerer fc de Barrosll2009h . 

We have extensively discussed the signatures of EMP 
and metal-free galaxies expected from the model. In current 
and future observational data, we can obtain a sample of 
galaxies at high- z by the standard drop-out technique. If we 
can select EMP and metal-free candidates from it by broad- 
band data, it is very useful to select the target for follow-up 
spectroscopy. The bluest galaxies in rest-frame UV are a 
good target (Figures 11 and 12). An example criterion in 
the colour with the HST/WFC3 is J 125 - Hieo < -0.15 
for z ~ 7-8. The blue gal axies with the UV spe ctral slope 
P = -3 ± 0.2 reported bv lBouwens et all l|2010ah are likely 
to be EMP or metal-free although we cannot deny the possi- 
bility that these galaxies have normal metallicity if they are 
very young (~ 1 Myr) and have a large LyC escape fraction 
(/esc > 0.5) (Figure 11). Rest-frame optical colours are also 
useful because there are some metallicity indicators such as 
the strong Balmer jump, Balmer series lines, and oxygen 
lines (Figures 13 and 14). For example, the criteria with the 



http: / /www. stsci.edu/jwst/instruments/nirspec/sensitivity/R100. 



JWST/NIRCAM colour are F277W-F444W > for z ~ 6 
and F277W-F444W < for z ~ 8. 

After selected the candidates, we will perform follow-up 
spectroscopy for them, and then, measure the metallicity. At 
this point, equivalent widths (EWs) of emission lines become 
useful as an indicator. We have examined EWs of Lya, He n 
A1640, Ha, H/3, and [O m] A5007 as a function of metallicity, 
star formation duration, and LyC escape fraction (Figures 
6-10). For the duration > 10 Myr and / csc = 0, the EMP 
and metal-free criteria of the rest-frame EW are as follows: 
EW(Lya) > 230 A, EW(He II) > 1 A, EW(Ha) > 1900 
A, EW(H/3) > 300 A, and EW(0 ill) < 20 A. Note that 
the Lya and He II criteria have a further uncertainty due 
to the resonant scattering in the ISM and IGM and due to 
the contribution of Wolf-Rayet stars, respectively. A crite- 
rion independent of the star formation duration and / eS c is 
the metal-to-hydrogen line ratio. The most easily detectable 
ratio is [O ill] A5007/H/3. If we find the ratio < 0.1 from a 
galaxy, it is EMP or metal-free. This ratio at z ~ 8 can be 
detectable by spectroscopy with the JWST/NIRSpec within 
a reasonable exposure time. In addition, we find that the 
Balmer jump can be also detected with the JWST/NIRSpec 
spectroscopy. 
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APPENDIX 

A machine-readable tables of emission line intensities rela- 
tive to H/3; Table Al is the no dust model and Table A2 is 
the dusty model. 



